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Abstract 

A calculation method based on an Interactive 
boundary-layer approach to multielement airfoils 
Is described and Is applied to three types of air- 
foil configurations with and without flap-wells In 
order to demonstrate the applicability of the 
method to general high-lift configurations. This 
method, well tested for single airfoils as a func- 
tion of shape, angle of attack and Reynolds number, 
Is here shown to apply equally well to two-element 
airfoils and their wakes, to a flap-well region, 
and to a three-element arrangement which Includes 
the effects of co-flowing regions, a flap well, and 
the wake of the elements. In addition to providing 
accurate representation of these flows, the method 
Is general so that Its extension to three- 
dimensional arrangements is likely to provide a 
practical, accurate and efficient tool to assist 
the design process . 

1 .0 Introduction 

The design of multielement airfoils for high 
lift requires consideration of a range of configu- 
rations so that care must be taken to ensure that 
the essential experiments and calculations can be 
undertaken with acceptable cost as well as accu- 
racy. * 2 We are concerned here with the devel- 
opment of a calculation method which meets this 
requirement and Is able to represent the flow over 
and between the Individual airfoils with consider- 
ation of flap wells and wakes. These requirements 
Imply the need for a method which has an economical 
and accurate numerical solution procedure; a flex- 
ible turbulence model to represent wall boundary 
layers, wall Jets and wakes, and the wake of the 
last element; and the ability to represent the 
separated flows associated with the upper surface 
and the flap well. In addition, the preferred pro- 
cedure should be readily extendible to deal with 
three-dimensional components such as wings and 
empennage. 

In recent years there has been a renewed Inter- 
est In experimental work on high-lift systems. 
Extensive measurements have been reported by 
Nakayama 0 , Alemdaroglu summarized by Nakayama 4 , and 
Valarezo et al. 5 The data of Nakayama are for a 
three-element airfoil with a leading-edge slat and 
for a single-segment flap; they were obtained at 
NASA Langley’s Low Turbulence Pressure Tunnel 
( LTPT) . Those of Alemdaroglu are essentially for 
the same but smaller model and were obtained at the 
low-speed wind tunnel of California State Univer- 
sity, Long Beach. The data of Valarezo et al. were 
also obtained at NASA Langley's LTPT and correspond 
to measurements at high Reynolds numbers. These 
data add to the previously obtained data on multi- 
element airfoils by van den Berg, 0 van den Berg 
and Oskam', Oskam et al. 0 , Omar et al, 9 »'° 
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and Olson and Orloff 11 which allow the validation 
of computer programs to analyze high-lift systems. 

Several airfoil-analysis and design algorithms 
have been developed In the last decade and have 
been based on one of two approaches: numerical 
solutions of the Reynolds-averaged Navler-Stokes 
equations or solutions of the Interaction between 
Invlscld and boundary-layer equations. The former 
approach Involves the numerical solution of ellip- 
tic equations so that Information travels In all 
directions through pressure, velocity and viscous 
and turbulent stress gradients. As a result, the 
solution method requires simultaneous processing 
of the pressure and velocity components and stress 
tensor throughout the flowfleld and this, In turn, 
Implies a trade-off between accuracy and cost which 
tends to limit the validity of this approach. This 
limitation Is a function of computers and program- 
ming methods, and these are likely to Improve with 
time so that solutions of the Navler-Stokes equa- 
tions, with proper consideration of momentum con- 
servation In two directions together with longitud- 
inal diffusion, are likely to be a major component 
of design methods of the future. The combination 
of the largest mainframe computers and unstructured 
and multigrid techniques has already been shown to 
be very powerful as, for example, by Mavrlplls 12 , 
Rogers et al.' 3 , and Barth.' 4 

The present approach Is based on the Interac- 
tive boundary-layer approach which has been tested 
extensively for single-element airfoils, as des- 
cribed for example In References 15 to 17. These 
papers have shown that this approach can represent 
accurately, and with low cost, the flows around a 
number of airfoil geometries, with angles of attack 
up to and beyond that of stall, and Including re- 
gions of separated flow which may cause transition 
from laminar flow. It Is based on solutions of 
Invlscld and boundary-layer equations with a sur- 
face and wake blowing velocity obtained from the 
Hilbert Integral and ensuring Interaction between 
the calculated Invlscld and viscous flows. It has 
also been shown that, In extended form. It Is able 
to represent the three-dimensional flows over 
wings 10 and, therefore, meets many of the require- 
ments for a design method as discussed above. 

It should be emphasized that alternative Inter- 
active methods have been reported, for example by 
Veldman 10 , LeBal leur 20 * 21 , Williams and Smith 22 , 
Drela and Giles 23 , and Kusunose, Wlgton and 
Meredith. 24 The last of these deserves further 
research In the context of this paper since It has 
been reported recently and has been applied to 
multielement airfoils. It uses a finite-element 
full-potential code to compute the outer flow with 
a modified streamline H-grld, and solves the 
boundary- layer equations In Integral form with 
modeling similar to that of Bradshaw and Ferris. 25 
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The results encompass single, two- and three- 
element airfoils at angles of attack up to around 
13 degrees and are In close agreement with mea- 
surements. It Is also evident that the method Is 
cost efficient and Is already part of a design 
method. It can be expected, however, that the 
turbulence model will be less successful where 
pressure gradients are severe and will not cope 
well with Inter-element flows where there Is a 
distinct velocity maximum. In addition, and per- 
haps of greatest Importance, attempts to extend 
Integral boundary-layer methods to three- 
dimensional flows have not been successful. 

This paper Is concerned with the extension of 
the Interactive boundary-layer method of Refs. 15 
to 17 to represent multielement airfoils where the 
flows between airfoils, flap wells and the possible 
influence of the overall wake In all elements are 
new features. The computational Investigation was 
carried out In three parts which are reflected In 
the presentation of results and correspond to two- 
element airfoils with emphasis on the flow between 
elements and the wake, a single-element airfoil 
with a flap well where the calculation of the flow 
In the flap well Is the major novelty, and the 
combination of these features In three-element 
airfoils which Involve a flap well. This results 
section Is preceded by descriptions of the Inter- 
active and solution procedures and followed by 
concluding remarks. 

2.0 Interactive Boundary-Layer Method 

The Interactive boundary-layer method makes use 
of the panel method of Hess and Smith 26 and a 
solution of the boundary-layer equations In which 
the turbulence model Is given by the algebraic eddy 
viscosity (e m ) formulation of Cebecl and Smith. 27 
With b denoting 1 + c m /v, the continuity and 
momentum equations can be written as 


&u e (x) 
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This Inverse boundary-layer formulation Is appro- 
priate to airfoils and to those parts of airfoils 
without surface discontinuities such as flap wells. 
Where flap wells occur, a different formulation of 
the Inverse procedure Is required, and the formula- 
tion used here Is described below. 


The calculation of the flow In the flap-well 
region Is similar to that over a backward facing 
step. A large portion of the flow separates 
Immediately after the sudden change of the geom- 
etry, and the size of the reversed-f low region 
depends mainly on the step height, on the gap, and 
the overhang. The flow reattaches and gradually 
recovers downstream In the flap-well region or In 
the wake. The calculation of flows of this kind 
Is difficult, and potential theory Is not adequate 
because of the singularity that occurs at the 
geometry discontinuity and the strong viscous 
effects In the separated flowfleld. Thus, an 
Initial distribution of displacement thickness Is 
assumed and the relaxation formula 

(4*) V+1 = («*) V [1 ♦ « - 1)] (6) 

U e1 

Is used In the Inverse method to replace the 
Hilbert Integral formulation of the external 
boundary condition. The new edge boundary condi- 
tions are given by Eq. (3b) and Eq. (6), where u ev 
and Up^ correspond to the external velocities com- 
puted by the boundary layer and Invlscld methods, 
respectively, and u Is a relaxation parameter. 
At the end of the flap-well region, the solution 
procedure reverts to the Hilbert-Integral approach. 

2.2 Turbulence Node! 


au ay 
ax ay 


0 


^ ♦ V & - u 
ax ay e 


du, 

dx 


* V 


3 /h 

ay (b ay* 


0 ) 

( 2 ) 


The turbulence model used to represent the flow 
on the airfoil may be expressed In terms of the 
Cebecl and Smith eddy-viscosity formulation, 
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In the absence of mass transfer, the boundary con- 
ditions for the above equations on the airfoil are: 


u = v = 0, 

y = 0 

(3a) 
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and In the wake, where a dividing line at y = 0 Is 
required to separate the upper and lower parts of 
the Invlscld flow and In the absence of the normal 
pressure gradient, they are: 

y -* ±co r u •* u e (x); y = 0 , v = 0 (4) 


2.1 Interaction law 

To perform the calculations for flows with 
separation. It Is necessary to use an Inverse 
procedure and compute the external velocity as 
part of the solution. Here we use the formulation 
discussed In Ref. 15 and write the edge boundary 
condition as 


u e (x) = ug(x) ♦ *u e (x) (5a) 

with 6u e (x) computed from the Hilbert Integral 
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The condition used to define y c Is the continuity 
of the eddy viscosity so that Eq. (7a) Is applied 
from the wall outward (Inner region) until Its 
value Is equal to that given for the outer region 
by Eq. (7b). The expression ytr represents the 
transition region and Is given by 

T = 1 - exp[G(x - x. ) J jp] (8b) 

x tr e 

Here xt r Is the location of the beginning of tran- 
sition and G Is defined by 
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where the transition Reynolds number R x . 
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The location of the onset transition Is 
obtained from Michel’s formula, given by Ref. 28, 


R e » 1.174 (1 + 22 xM) rO-46 


(9) 


with Vmln denoting the location where the velocity 
Is a minimum. 

The eddy-viscosity model for the flap-well, 
like all expressions for turbulent flows, Is empir- 
ical and was first tested for flow over a backward 
facing step before Its application to the present 
problem. While the agreement with backward facing 
data was satisfactory. It should be tested further 
and possibly replaced by "better" expressions or 
models. 

2.3 Solution Procedure 


When flow separation takes place upstream of the 
transition location predicted by this formula, 
transition Is assumed to coincide with the location 
of separation. 


In the flap-well region, the above formulas are 
modified so that 
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Here c(^) denotes the eddy viscosity corresponding 
to the velocity profile above the separated region 
and c^F) includes the total region from the wall. 
The expressions for c ft, 1 ) and cfoF' are given by: 
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where y 0 Is the location of u = 0, x Is a relaxa- 
tion parameter (usually around 10), L is a charac- 
teristic length, and x 0 Is the beginning of the 
flap-well . 


In the wake the corresponding expressions are: 


In general. It Is convenient to solve the equa- 
tions of the previous section In transformed vari- 
ables. The Falkner-$kan transformation defined by 

u e 1/2 1/2 

n = (^) y, = (u e »x) r/<? f(x,n) (15) 

Is used here and, with the usual definition of 

stream function, Eqs. (1) - (3) lead to 

(bf") 1 * 1 (m * 7)ff " + m(1 - f ,2 ) 
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n * 0, f * f’ « 0, t> = n e » f* » 1, (17) 

where primes denote differentiation with respect to 
n and 


A slightly modified form of this transformat Ion 
Is used when the calculations are performed In the 
Inverse mode by replacing u e with a reference 
velocity u 0 , that Is, 


Y = /u Q /vX y, * = /U Q vx F(x,Y) (18) 


so that the continuity and momentum equations and 
their boundary conditions, given by Eqs. (16) and 
(17) become, 
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The boundary conditions corresponding to Eq. 
(5b) are obtained by applying a discretization 
approximation to the Hilbert Integral, Eq. (5b), 


where UmH.e. the eddy viscosity at the trail- 
ing edge computed from Eqs. (7) and (10) and (c m ) w 
Is the eddy viscosity In the far wake given by the 
larger of 
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where the subscript 1 denotes the x-statlon where 
the Inverse calculations are to be performed, c^j 
Is a matrix of Interaction coefficients, and 0 is 
given by 0 = u e 6*. Further details are available 
In Ref. 15. In terms of transformed variables, the 
parameter D can be written as 


o - <g) 1/2 (Y e w - F ) (22) 
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and the relation between the external velocity w 
and displacement thickness 6* provided by the 
Hilbert Integral can then be written as 

Y . Y e , w * (^) 1/Z (Y e - F) ♦ g, (23) 

where 
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In the flap-well region, Eq. (23) Is replaced by 

u Q 1/2 

Y.Y e . F e -w[Y e -(^) **] (25) 


The corresponding boundary conditions In the wake 
are 


Y * Y_ e# F* - w; Y =* 0, F - 0; 

Y « Y e , F' » w, 

with w now given by 

w = c 11 (w(Y e - Y_ e ) - (F e - F e )] (^) 1/2 ♦ g r 

0 (26) 

The solution of the above transformed equations Is 
obtained by Keller's box method, as described In 
Ref. 15. Where separation occurred, the convective 
term u(au/ax) was set to zero and this assumption 
proved to be satisfactory for the flow on the air- 
foil. The larger regions of separation associated 
with the flap well and the near wake required an 
additional Iteration scheme based on a continuation 
method. Since a linearized form of the boundary- 
layer equations Is being solved. It Is necessary 
that the calculations at station have Initial 
profiles which are usually assumed to correspond 
to those at a previous x-statlon, x^. With 
Increasing flow separation, the effect of the Ini- 
tial profiles on the solutions at Increases 

and can lead to breakdown. A remedy to this prob- 
lem Is to define the velocity profile at x^ to 
be of the form 

u = u re f + n(u a - u re f ) (27) 

where u a denotes the velocity profile at x^i and 
u re f to a profile which allows solutions to be 
obtained at xy The Iteration process at x 1 pro- 
ceeds with values of n ranging from 0 to unity. 


The sequence of the calculations Is as follows. 
The panel method provides an external velocity 
distribution based on a body shape In which the 
flap-well region Is assumed to be absent. The 
Interactive boundary-layer approach leads to solu- 
tions on the upper surface from the stagnation 
point through the regions of laminar, transitional 
and turbulent flow to the trailing edge. Simi- 
larly, it provides results for the lower surface 
up to the beginning of the flap-well. A displace- 
ment thickness distribution, 4*(x), Is assumed 
In the flap well and, with the continuation method 
described above and with the Initial velocity pro- 
file similar to that of a backward- facing step, 
calculations proceed to the trailing edge. With 
the upper and lower surface velocity profiles at 
the trailing edge, the calculations are extended 


Into the wake. As a consequence of the above, a 
blowing velocity Is available on the airfoil and 
In the wake. In the flap-well region, the blowing 
velocity v n Is defined by 


v * ~ (u 6*) 
n dx ' e m' 


(28) 


where 6^ =* A* - if Here 5 t corresponds to the 
body shape assumed to exist over the flap-well. 
Elsewhere, the blowing velocity Is given by 
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With the blowing velocity distribution known 
everywhere In the flowfleld, a new distribution of 
external velocity Up^(x) Is obtained from the 

panel method. As before, the boundary- layer jfcolur 
tlons on the upper and lower surfaces of the air- 
foil are obtained with the Hilbert Integral In 

which the edge boundary condition, Eq. (21), Is 
now written as 
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with k Indicating the Iteration cycle. Equation 

(23) also changes to 
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In the flap-well, with u ev known from the 
previous flap-well calculation, a new ^-distri- 
bution is available from Eq. (6) and Is used to 
obtain solutions up to the trailing edge. This 
sequence of calculations Is repeated for the whole 
flowfleld until convergence Is achieved. The con- 
tinuation method discussed before Is Involved 
within this sequence of calculations where 
necessary. 


3.0 Results and Discussion 
3.1 Two-Element Airfoil 

In a previous study, a similar Interactive 
approach was used to compute the performance char 
acterlstlcs of three two-element airfoils. 29 
The 1 nvl sc Id flow solutions were obtained By the 
conformal-mapping method of Halsey 30 , rather 
than the panel method used here, and the viscous- 
flow calculations were performed without ac cou nting 
for the wake effects, either behind the main air- 
foil or the flap. The calculation method provided 
results which agreed well with experimental Infor- 
mation within the accuracy of the measurements up 
to an angle of attack which was sufficiently small 
so that there was either no or very small separa 
tlon on the airfoil and the gaps between the ele- 
ments were comparatively large. In this way, the 
difficulties In computing the wake of each airfoil 
and accounting for the merging of the shear layers 
between the airfoil and the flap, and extending the 
range of the computational method to higher angles 
of attack were postponed to a later time. 
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In the studies reported here, we first per- 
formed calculations with the present method which 
did not Include the wake effect and compared the 
results with those obtained with the earlier 
code 29 with Its different Invlscld flow method. 
After ensuring that the results of both codes were 
essentially the same, the wake effects were Intro- 
duced Into the present method and calculations 
were repeated for the three two-element airfoils 
to Investigate the role of the wake effect on the 
solutions . 

The first results of this paper are for two- 
element airfoils for which corresponding experi- 
mental Investigations have been reported by Van 
den Berg 8 , by Omar et al. 9 »^° and Olson and 
Orloff. 1 ' There are no flap wells In these 
arrangements and the novel features of the calcu- 
lations are the flow between elements and the 
Inclusion of the wakes. 

Figures 1-3 present the results of the data of 
Van den Berg 8 , also discussed by Van den Berg 
and Oskam 7 and Oskam et al. 8 , which correspond 
to a supercritical main airfoil (NLR 730) with a 
flap of 32% of the main chord at a deflection angle 
of 20 degrees. Measurements of surface pressure 
and velocity profiles were obtained a chord 
Reynolds number of 2.51 x 10 5 and for angles of 
attack of 6 and 13.1 degrees, the latter recognized 
as the highest angle which corresponded to fully 
attached flow. Lift coefficients were deduced for 




(b) 

Fig. 1. NLR 7301 wing with flap. Calculated and 
measured: pressure distribution at (a) a = 6°, 

(b) a = 13. r. 



Fig. 2. NLR 7301 wing with flap. Calculated and 
measured lift coefficients. 



(a) 



Fig. 3. NLR 7301 wing with flap. Calculated and 
measured: local skin-friction coefficients, Cf, 

and momentum thicknesses, e/c, on the upper wing 
surfaces at (a) a = 6°, (b) a * 13.1°. 

five angles of attack. The airfoil arrangement Is 
shown on Fig. la and It Is evident that the dis- 
tributions of pressure coefficients are In close 
agreement with the measured values over the sur- 
faces of both elements. The results at the angle 
of attack of 13.1 degrees confirm the absence of 
separation and the lift coefficients of Fig. 2 
that stall occur at an angle larger than this 
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value. The calculated variation of lift coeffici- 
ent with angle of attack Is close to the measure- 
ment and some 3% lower than that calculated without 
consideration of the wake. It Is to be expected 
that this difference will Increase with angle of 
attack and particularly as separation occurs and 
expands over the upper tralllng-edge region, and 
this trend Is evident In the figure. 

figures 3a and 3b show the variations of 
momentum thickness and skin-friction coefficient 
with chord distance over the main airfoil. The 
agreement between calculated and measured results 
Is remarkable for both angles of attack, the only 
significant discrepancies being In the skin- 
friction coefficient In the upstream region of max- 
imum rate of change. Again, the results at 13.1 
degrees confirm the absence of separation, 

although this result conveys little about the flow 
on the second element. 

Figure 4 presents the results for a NASA super- 
critical airfoil, 24 M In length, with a 7" flap at 
a deflection angle of 20 degrees. The experiments 
were carried out In the 36 x 96 In. wind tunnel of 
the Boeing Research Laboratories at a Mach number 
of 0.2 and have been documented by Omar et al. • 
The pressure-coefficient distributions of Figs. 4a 
and 4b correspond to angles of attack of zero and 



(a) 



Fig. 4. NASA supercritical two-element airfoil. 
Calculated and measured: (a) pressure distribution 
at a = 0°, (b) pressure distribution at a = 8.93°, 
The results of Ref. 29 without wake effect. 


8.93 degrees with the measured and calculated val- 
ues again within experimental uncertainty. The 
wake had no effect at zero angle of attack, as 
expected, and had a slight effect on the pressure 
coefficients at the 8.93 degree angle, although the 
Improvement on the main airfoil Is coupled with an 
apparent lack of Improvement on the flap. These 
results are reflected In the lift coefficients of 
Fig. 5 where the calculated results with the wake 
effect are In better agreement with data than those 
without the wake effect. 



results of Ref. 29 are without wake effect. 

The third two-element airfoil corresponds to 
that Investigated by Olson and Orloff” which 
Involves a NACA 4412 airfoil with a chord length 
of 0.9m upstream of a flap which has the section 
of a NACA 4415 airfoil and a chord of 0.36m with a 
deflection angle of 10°. Figure 6 shows the mea- 
sured and calculated surface pressure distributions 
for a Reynolds number of 1.3 x 10 b at an angle 
of attack of 2.2°. As In the first case, there Is 
a slight Improvement over the results obtained 
without the wake which, again, may be due to the 
absence of flow separation. 



It Is evident from the comparisons of Figs. 1-6 
that the flow between the airfoil elements and the 
wake have been satisfactorily Incorporated In the 
Interactive procedure with results which are virtu- 
ally Identical to the measurements. It should, 
however, be noted that the effect of merging of the 
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boundary layers In these three configurations Is 
negligible. 

3.2 Single Airfoil with a Flap-Well 

An appraisal of the Interactive boundary-layer 
procedure, as applied to the flow In and around a 
flap well, required corresponding measurements and 
a parallel experimental program which was carried 
out at the California State University, Long Beach, 
with the single airfoil arrangement of Fig. 7. The 
experiment described In Ref. 4 used this 12-Inch 
chord airfoil, which corresponds to the main ele- 
ment of the three-element configuration tested In 
the low-turbulence tunnel of the NASA Langley 
Research Center, and described by Nakayama. 3 The 
chord Reynolds number was 0.5 x 10 6 , and surface 
pressures were measured for angles of attack up to 
14 degrees with local velocity Information In the 
flap well at an angle of attack of 5 degrees. Fig- 
ure 8 presents the measured and calculated surface- 
pressure distributions for angles of attack of 5, 
8 and 12 degrees with transition tripped at 0.25c, 
and the agreement Is generally good. The calcu- 
lated upper-surface pressure peak close to the 
leading edge reflects the better spatial resolution 
of the Interactive method, and close to the trail- 
ing edge there are some small disagreements which 
may stem from the flap-well results. Neverthe- 
less, the pressure-coefficient distributions 
represent closely the measurements in the flap- 
well with the near constant values Indicating the 
region of recirculation followed, as can be seen, 
by a rapid Increase In pressure coefficient after 
reattachment. 



Fig. 7. Single airfoil with flap-well cut. 


Figure 9 shows calculated profiles of stream- 
wise velocity, streamlines and the distribution of 
skin-friction coefficients within the flap well. 
The profiles differ Increasingly from the measure- 
ments In the near-wall region, although reattach- 
ment occurs at about the same location. The 
reasons for the discrepancy are likely to be 
associated with the turbulence model that Is used 
here to represent a near-wall flow, which undoubt- 
edly Involves low- Reynolds-number characteristics 
together with a switch from a wall Jet In the nega- 
tive to a boundary layer In the positive direction. 
It Is unlikely that these local differences will 
affect greatly the outer region, particularly since 
the negative velocities and momentum are small, so 
that the calculated distribution of displacement 
thickness of Fig. 9b Is likely to be correct. The 



(a) 



Fig. 8. Comparison of calculated pressure distri- 
butions with data for R c =* 0.5 x 10 & at (a) 
a = 5.0, (b) a = 0.0, (C) a = 12.0. 

distribution of skin-friction coefficient, Fig. 9c, 
Indicates negative values up to reattachment at 
0.87c which Is In agreement with the velocity 
profiles . 

3.3 Three-Element Airfoils 

With the positive results of Figs. 1 to 9 for 
three two-element airfoils and for an airfoil with 
a flap-well, It Is appropriate to consider the 
application of the Interactive boundary- layer 
method to a three-element airfoil with a flap well. 
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Fig. 9. Results for R c - 0.5 x 10 & at a = 5.0. 
(a) Velocity profiles In the flap-well region, (b) 
Recirculation streamline and the location of the 
displacement thickness, (c) Calculated local skin- 
friction coefficient In the flap-well region. 


The chosen configuration Is shown on Fig. 10 and 
corresponds to the high-lift model, tested In the 
NASA, Langley and California State wind tunnels at 
Reynolds numbers of 5 x 10 6 and 0.5 x 10°, respec- 
tively. 3 * 4 The slat deflection angle was -30 
degrees and the flap deflection angles 15 and 30 
degrees with angles of attack of 4 to 20 degrees. 
The measurements were made by a combination of hot- 
wire and laser-vel oclmetry techniques, the latter 
was primarily used In regions of separated flow. 

I 

I 
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Fig. 10. The three-element airfoil with analytical 
and experimental fairing, i 5 =» -30°, if » 15 p . 

Calculations were Initially made on smooth 
bodies without explicitly considering the flow In 
the flap-well region. Also, because the potential 
flow theory predicts flow singularities at the 
discontinuity of the airfoil geometry, the sharp 
corner of the slat and the flap-well cut out of 
the main airfoil were smoothed to prevent solu- 
tions from breaking. Figure 10 shows the modified 
geometry of this airfoil with the flap-well fairing 
and the rounded slat used In calculations. The 
so-called "experimental fairing" refers to the 
dividing streamline which was determined from mea- 
surements, while the "analytical fairing" was drawn 
arbitrarily. Figure 11 shows the velocity vectors 
for a particular combination of gap and overhang 
and the position of the streamline dividing the 
recirculating flow behind the flap-well step from 
the outside flow, as determined from the mean- 
velocity vector data. The position of this divid- 
ing streamline Is Important since It corresponds 
roughly to the equivalent smooth body with pressure 
distribution close to the real one. Hence, this 
dividing streamline was used as the "experimental 
fairing" in the calculations. Figure 12 shows the 
surface-pressure distributions on the slat, main 
airfoil, and flap at three angles of attack (4°, 
12° and 16°) for the configuration with the experi- 
mental fairing, and Fig. 13 shows the corresponding 
distributions at the same angles of attack for the 
configuration with analytical fairing. Overall, 
the calculated results agree well with experi- 
mental data except for the slat at low angles of 
attack and the pressure peak on the main airfoil. 
This discrepancy may be caused by the differences 
between the assumed slat shape and the real one. 
In the flap-well region, the results are In better 
agreement with measurements when the fairing Is 
close to a real streamline. 

Finally, the method of Section 4.2 was applied 
to Include the calculation Inside the flap-well. 
The results, shown In Fig. 14, agree well with 
measurements for all the cases Indicated above and 
Including angles of attack up to 20°. This con- 
firms that It Is unnecessary to make a priori 



Fig. 11. Velocity vectors In the flap-well region 
by laser velocimeter. 
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(C) 

Fig. 12. Pressure distribution on the three- 
element airfoil with experimental fairing for 
4 f = 15° and R c = 5 x 10 6 , (a) o = r, (b) a = 12°, 
(C) a = 16°. 

assumptions about the fairing shape, and allows 
for further detailed Investigations of the recirc- 
ulation flow In the flap-well, such as the gap and 
overhang effects. Figure 15 shows the variation 
of lift coefficient with angle of attack, confirm- 
ing that the present calculation method leads to 
values which are In close agreement with 
experiment . 


Fig. 13. Pressure distribution on the three- 
element airfoil with analytical fairing for if 
= 15° and R c = 5 x 10 6 , (a) a * 4°, (b) a 
. 12°, (c) a * 16°. 

Comparison of pressure coefficients for a flap 
deflection angle of 30°, Fig. 16, allows similar 
conclusions to be drawn to those of the previous 
paragraph. Also, the calculated lift coefficients 
shown In Fig. 17 are very close to measurements. 

4.0 Concluding Remarks 

The results and discussion of the previous 
section show that the present Interactive method, 
with Its consideration of the flap-well region and 



















the wakes, leads to pressure coefficient distribu- 
tions and values of lift which are In good agree- 
ment with experiment for single and multielement 
configurations with angles of attack up to 20°, 

The present method neglects the confluent 
boundary layers, which become Important as the 
distance between the main airfoil and the flap 
becomes smaller or where the shear layer thickens 
with angle of attack. Also, close examination of 
the turbulence model In the vicinity of this merg- 
ing region Is desirable. The Incorporat Ion of 
Improvements to consider these aspects Is likely 
to lead to a more generally applicable method. 



Fig. 17. Variation of the lift coefficient with 

angle of attack for the three-element airfoil with 

6f = 30* and R c * 5 x 10 6 . 
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